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Abstract 

A semicoherent system can be described by its structure function or, equivalently, 
by a lattice polynomial function expressing the system lifetime in terms of the compo- 
nent lifetimes. In this paper we point out the parallelism between the two descriptions 
and use the natural connection of lattice polynomial functions and relevant random 
events to collect exact formulas for the system reliability. We also discuss the equiv- 
alence between calculating the reliability of semicoherent systems and calculating the 
distribution function of a lattice polynomial function of random variables. 

Keywords: Reliability, semicoherent system, lattice polynomial function, weighted lattice 
polynomial function. 



1 Introduction 

Consider a semicoherent system made up of nonrepairable components. Such a system can 
be described by its structure function, which expresses at any time the state of the system in 
terms of the states of its components. Equivalently, the system can be described by a lattice 
polynomial (l.p.) function which expresses the system lifetime in terms of the component 
lifetimes. 

In this paper, we point out the formal parallelism between both descriptions, we collect 
exact formulas for the system reliability, and we show that calculating the reliability of 
semicoherent systems is equivalent to calculating the distribution function of an l.p. function 
of random variables. 

* Corresponding author. 
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We also consider the more general case where there are collective upper bounds on life- 
times of certain subsets of components imposed by external conditions (such as physical 
properties of the assembly) or even collective lower bounds imposed for instance by back- 
up blocks with constant lifetimes. In terms of lifetimes, such systems can be described by 
weighted lattice polynomial (w.l.p.) functions. In terms of state variables, we will see that 
a "weighted version" of the structure functions is required. 

This paper is organized as follows. In §2 we discuss the parallelism between the descrip- 
tion of semicoherent systems by structure functions and by the corresponding l.p. functions. 
In particular, in §2.3^ Theorem [2] uses the natural connection between lattice polynomial 
functions and relevant random events to establish a centrally important relation between 
the lifetimes of the system and its components. In §3 we yield exact formulas for the system 
reliability in case of independent arguments and in general. In turn, those formulas make 
it possible to provide exact formulation of reliability parameters such as the mean time-to- 
failure of the system. In §4 we generalize our results by considering lower and upper bounds 
on lifetimes of certain components. Finally, in §5 we examine how our results can supply 
exact formulas for the distribution and moments of w.l.p. functions of random variables. 

For any numbers a, /3 G M := [—oo, oo] and any subset A C [n] := {1, . . . , n}, let e^'^ 
denote the characteristic vector of A in that is, the n-tuple whose ith coordinate is 

(3, if i E A, and a, otherwise. Also, the Li-norm of any binary vector x G {0, 1}" is denoted 



2 Structure function and l.p. function 

In this section we recall the main concepts and results related to structure functions of semi- 
coherent systems. We also point out the parallelism between the description of a system by 
its structure function and the description of this system by an l.p. function of the component 
lifetimes. 



2.1 Structure function 

Consider a system consisting of n components that are interconnected. The state of a 
component i G [n] can be represented by a Boolean variable Xi defined as 




1, if component i is functioning, 
0, if component i is in a failed state. 



For simplicity, we also introduce the state vector x 

The state of the system is described from the component states through a Boolean func- 
tion : {0, 1}" — i> {0, 1}, called the structure function of the system and defined as 




1, if the system is functioning, 
0, if the system is in a failed state. 



We shall assume throughout that the structure function (p is nondecreasing (the system is 
then said to be semicoherent) and nonconstant, this latter condition ensuring that 0(0) = 
and 0(1) = 1, where := (0, . . . , 0) and 1 := (!,...,!). For a background on semicoherent 
systems, see for instance the monographs by Ramamurthy |14j and Rausand and H0yland 



As a Boolean function, the structure function (j) can also be regarded as a set function 
V : 2'"] — s> {0, 1}. The correspondence is straightforward: We have v{A) = (/)(e^^) for all 
A C [n] and 

0(x)= viA)Ylx, n (1-^^)- (1) 

AC[n] i&A i&[n]\A 

We shall henceforth make this identification and often write 0„(x) instead of (/'(x). Clearly, 
the structure function is nondecreasing and nonconstant if and only if its underlying set 
function v is nondecreasing and nonconstant. 

We also observe that, being a Boolean function, the function 0^ has a unique expression 
as a multilinear polynomial in n variables, 

(t>M) = J2 ^viA) n X, (2) 

AC[n] i&A 

(see for instance Hammer and Rudeanu [7]), where the set function niv ■ 2'"! — > Z is the 
Mohius transform of f , defined by 

mM)= E(-l)'^'"""^(^)- 

Another concept that we shall often use in this paper is the dual of the set function f , 
that is, the set function v* : 2^ {0, 1} defined by v*{A) = 1 - v{[n] \ A). 

By extending formally the structure function 0^ to [0, 1]" by linear interpolation, we 
define the multilinear extension of 4>y (a concept introduced in game theory by Owen [13]), 
that is, the multilinear polynomial function 0^ : [0, l]** [0, 1] defined by 

?.(x)= E viA)l[x, n (1-^0- (3) 

AQn] ieA ie[n]\A 

Now, by combining the concepts of Mobius transform, dual set function, and even the 
"coproduct" operation 11, defined by IljXi = 1 — nj(l — Xj), we can easily derive various useful 
forms of the structure function. Each of these forms is a polynomial expression of the function 
(j)^ and hence, when formally regarded as a function from [0, 1]" to [0, 1], it identifies with 
the corresponding multilinear extension 0^; see also Grabisch et al. [5J. Table [1] summarizes 
the best known forms of the structure function and its multilinear extension. 

2.2 L.p. function 

For any event E, let Ind(i?) represent the indicator random variable that gives 1 if i? occurs 
and otherwise. For any z G [n], we denote by Tj the random time-to-failure of component 
i and we denote by Xi{t) = Ind(Tj > t) the random state at time t ^ of component i. For 
simplicity, we introduce the random time-to-failure vector T = (Ti, . . . , r„) and the random 
state vector X(t) = {Xi{t), . . . at time t ^ 0. We also denote by Ts the random 

time-to- failure of the system and by Xs{t) = Ind(T5' > t) the random state at time t ^ of 
the system. 

The structure function 0„ clearly induces a functional relationship between the variables 
Ti, . . . ,Tn and the variable T5. As we will see in Theorem |51 Ts is always an l.p. function 
of the variables Ti, . . . ,T„. Just as for the structure function, this l.p. function provides a 
complete description of the structure of the system. 
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Name 


0^(x) and 0„(x) 


Pri m a,l form 


V v(A) Y\ X TT (1- x ) 

AC[n] ieA ie[n]\A 


T)llJ^l TOTTTl 
A—y Lldl IVJI 111 


1 - V v*( A\ TT T- TTri - r \ 
^ 2^ ^ \-^) 11 -^-^ 111-"- •^■^1 

AC[n] iG[n]\A i&A 


PnTTipl A/Tonmc; ToyTTi 

J. 1 lllldl IVivJUlLlo IWl 111 


m { A\ TT T- 
/ , iit-v\^) 11 Xi 

AC[n] i€A 


T)npil A/Tonint; tottti 

J--' Lldl IVl-WUlLlO IWl 111 


\ rn . f /l"! TT T- 
AC[n] ieA 


Disjunctive normal form 


II v{A)\{x, 

Ac[ji] ieA 


Conjunctive normal form 


n ^*(^)ii^. 

AC[n] ieA 



Table 1: Various forms of the structure function and its multilinear extension 

Let us first recall the concept of l.p. function of real variables; see for instance Birkhoff [U 
§11.5] and Gratzer P, §1.4]. Let L C M denote a totally ordered bounded lattice whose lattice 
operations A and V are respectively the minimum and maximum operations. Denote also 
by a and h the bottom and top elements of L. 

Definition 1. The class of lattice polynomial (l.p.) functions from L" to L is defined as 
follows: 

(i) For any k E [n], the projection t i— > is an l.p. function from to L. 

(ii) If p and q are l.p. functions from to L, then p /\ q and pM q are l.p. functions from 
L" to L. 

(iii) Every l.p. function from to L is constructed by finitely many applications of the 
rules (i) and (ii). 

Clearly, any l.p. function p : L" — * L is nondecreasing and nonconstant. Furthermore, 
it was proved (see for instance Birkhoff p^, §11.5]) that such a function can be expressed 
in disjunctive and conjunctive normal forms, that is, there always exist nonconstant set 
functions w''' : 2["1 — > {a, h] and : 2^"^^ {a, b}, with w'^{0) = a and w'^{0) = b, such that 

p(t)= V A^«= A V^- (4) 

AC[n] ieA AC[n] iGA 

w'i{A)=b «)=(A)=a 

Clearly, the set functions w'^ and w'^ that disjunctively and conjunctively define the polyno- 
mial function p{t) in (j4]) are not unique. However, it can be shown [8] that, from among 
all the possible set functions that disjunctively define p{t), only one is nondecreasing. Sim- 
ilarly, from among all the possible set functions that conjunctively define p{t), only one is 
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nonincreasing. These special set functions are given by 



w\A)=p{ef) and = ^(6^4^. 

The l.p. function disjunctively defined by a given nondecreasing set function w : 2'"] — > {a, b} 
will henceforth be denoted Pw We then have 

pM= V A'^^= A V^- 

AC[n] i&A AC[n] i&A 

w{A)=b w*{A)=b 

where w* is the dual of w, defined as 

to* = 7 o V*, (5) 

the function 7 : {0,1} {«, being a simple transformation defined by 7(0) = a and 
7(1) = b. 

2.3 System descriptions 

The following theorem points out the one-to-one correspondence between the structure func- 
tion and the l.p. function that expresses Ts in terms of the variables Ti, . . . , r„. As lifetimes 
are [0, oo]-valued, we shall henceforth assume without loss of generahty that L = [0, 00], that 
is, a = and b = 00. We also make use of the transformation 7 as defined in (l5l). 

Theorem 2. Consider a system whose structure function (py : {0, 1}" — ^ {0, 1} is nonde- 
creasing and nonconstant. Then we have 

Ts = PwiTi,. . . ,Tn), (6) 

where w = 'y o v . Conversely, any system fulfilling ^ for some l.p. function p^, : ^ L 
has the nondecreasing and nonconstant structure function (p^, where v = o w . 

Proof. The proof mainly lies on the distributive property of the indicator function Ind(-) 
with respect to disjunction and conjunction, namely 

lnd{E V E') = lnd{E) V Ind(£;') and Ind(E A E') = lnd{E) A Ind(£;') 

for any events E and E'. Thus, for any t ^ we have 



Ind(p^(T)>t) = lnd( V A^*>^) 

V ACM ieA J 



v{A)=l 

= V Aind(r,>t) = II Y[x,{t) 

AC[n] i&A AC[n] ieA 

v{A)=l v{A)=l 

Hence, we have Ts = Pw{T) if and only if Xs{t) = 0t,(X(t)) for all t ^ 0, which completes 
the proof. □ 
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Remark 3. Since (py is a Boolean function, we can always replace in its expression each 
product n and coproduct 11 with the minimum A and the maximum V, respectively. Thus, 
Theorem [2] essentially states that (py is also an l.p. function that has just the same max-min 
form as but applied to binary arguments. More precisely, (/)y is similar to in the sense 
that 7 o 0„ = o (7, . . . ,7). 

We observe that many properties of the structure functions can be derived straightfor- 
wardly from the properties of the corresponding l.p. functions. Let us examine some of them 
(see for instance Rausand and H0yland [15, §3.11]): 

1. Boundary conditions. From the idempotency of p (that is, p{t, . . . , t) = t for all 
t e L), we immediately retrieve the idempotency of 0, that is, the boundary conditions 
0(0) = and 0(1) = 1. 

2. Internality. The internality property of p, namely 

n n 

1=1 1=1 

corresponds to the following internality property of 0: 

n n 
Y[Xi ^ 0(x) ^ Y[Xi. 



i=l i=l 

Note that, in both cases, internality results immediately from increasing monotonicity 
and idempotency. For instance, we have 

n / n n \ / n n \ n 

/\u = p( /\u,..., /\tA ^p{t) ^p(y u,...,y u) = yu. 

1=1 ^1=1 i=l ' ^1=1 j=l ' 1=1 

3. Pivotal decomposition. Consider the following median-based decomposition for- 
mula [H], which holds for any l.p. function: 



p{t) = median (^p(ai,t), U, p{bi,t)j, (7) 

where the ternary median function is defined as 

median(a;i, 3:2, X3) := (xi A X2) V (x2 A X3) V (^3 A Xi), 

and where (ai,t) (resp. (6i,t)) represents the vector t whose ith coordinate has been 
replaced with a (resp. b). From this formula we derive the following property of the 
structure function: 



0(x) = median(^0(Oj, x), Xj, 0(lj, X 

= 0(O„x) V (x, A0(l„x)) = 0(Oi,x)n(xin0(li,x 

= 0(Oi, X) + Xi0(li, X) - Xi0(Oi, X)0(li, X) 

and hence we retrieve the pivotal decomposition of the structure function, namely 

0(x) = Xi0(li,x) + (1 - a;i)0(Oi,x). 
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4. Structures represented by paths and cuts. From any nonconstant and nonde- 
creasing set function w : 2'"] define the set function u^, ■ 2'"] {a, ^} 
as 

, ib, if w(A) = b and w(B) = a hi all B Q A, 
Uw{A) = \ . (8) 

la, otherwise. 

The disjunctive and conjunctive representations of the l.p. function p^j having a mini- 
mal number of terms write (see Marichal [H Proposition 8]) 

pM= V A^^= A V^- (9) 

AC[n] ieA AC[n] ieA 

u^{A)=b u^*{A)=b 

Let us show that these representations are in one-to-one correspondence with the rep- 
resentations of the structure function by minimal paths and cuts. Recall that a path 
set P C [n] is a set of components which by functioning ensures that the system is 
functioning. Similarly, a cut set K C [n] is a set of components which by failing causes 
the system to fail. In other terms, P C [n] is a path set if v{P) = 1 and K C [n] is a 
cut set if v{[n] \ K) = 0. A path (resp. cut) set is minimal if it does not contain any 
proper path (resp. cut) set. 

It is known that if Pi , . . . , are the minimal path sets and Ki, . . . , Kg are the minimal 
cut sets, then 

r s 

M^) = II n ^» = n ]i 

j = l i(zp. j = l ifzKj 

The corresponding formulas for the l.p. function write 

r s 

pM = V a = a V 

j=i iePj j=i i&Kj 

and are exactly the "minimal" representations (jH]) oi Pw 

5. Extra component connected in series or parallel. Any l.p. function p : L"" ^ L 
fulfills trivially the following functional equations 

p{u Ati, . . . ,U Atn) = U Apiti, . . . ,tn), 

p{uy ti,. . . ,uy tn) = uy p{ti,. . . ,tn), 

for arbitrary u & L. These equations mean that connecting in series (resp. in parallel) 
any extra component to the system amounts to connecting that component in series 
(resp. in parallel) to each component of the system. The corresponding equations for 
the structure function are clear. We have 

(j){yxi,...,yxn) = 2/0(a;i, . . . ,x„), 
(t){yUxi,. . . ,yUxn) = yIl(f){xi,...,Xn), 

for arbitrary y E {0,1}. 



6. Dual structure. Recall that the dual structure function of a structure function (py 
is defined as 4>^{'x) = 1 — 0^(1 — x). From this definition, we derive immediately 
(j)^ = (j)y*, and hence from ([T]) we immediately retrieve the dual form of 0„ (i.e., the 
second expression in Tabled]). Using the dual set function w* of w, as defined in ([5]), 
we see that the corresponding l.p. function is the dual of p^j, namely p^ = pw*. 
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3 Exact reliability calculation 



The reliability function of component i is defined, for any t ^ 0, by 

R,{t) = Pr(T, >t) = Pr(X,(t) = 1) = E[X.,{t)], 

that is, the probabihty that component i does not fail in the time interval [0,t]. Similarly, 
for any t ^ 0, the system reliability function is 

Rs{t) = Pt{Ts >t) = PT{Xs{t) = 1) = E[Xs{t)], 

that is, the probability that the system does not fail in the time interval [0, t]. 

The mean time-to-failure of component i is defined as MTTFj = E[Tj] and similarly the 
mean time-to-failure of the system is defined as MTTF5 = E[r5]. These expected values can 
be calculated by the following formulas (see for instance Rausand and H0yland [151 §2.6]) 

MnFi= f\i{t)dt and MTTFs = [\sit) dt. 

J a J a 

In this section we yield the main known formulas for the system reliability function in 
the general case of dependent failures and in the special case of independent failures. We 
also provide some additional useful formulas. 



3.1 Dependent failures 

Dukhovny [2] found simple and concise formulas for the system reliability function in case 
of generally dependent variables Ti, . . . ,T„. We present them in the following theorem and 
we provide a shorter proof. 

Theorem 4. We have 

Rs{t) = E v{A)FT{X{t) = ef), (10) 



AC\n] 



Rsit) = l-Y. v*iA) Pr(X(t) = el-^J. (11) 

AC[n] 

Proof. By ([T]), we have 

Rs{t) = E[Mm)] = Y.<A)^\EMt) Ei^-Mt))] (12) 

AC[n] '-ieA ie[n]\A 

= <A)FTiXit) = ef), 

AC[n] 

which proves (|TOl) . Formula f|TT]) can be proved similarly by using the dual form of <py (i.e., 
the second expression in Table [I]) . □ 

Consider the joint distribution function and the joint survival function, defined respec- 
tively as 

F(t) = Pr(7^, ^t,\/ie [n]) and R{t) = Fi{Ti > t,\/i E [n]). 

By using the same argument as in the proof of Theorem HJ we obtain two further equivalent 
expressions of Rs{t). 
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Theorem 5. We have 



Rs{t) = mM)R{e7) 

AC[n] 

Rs{t) = 



(13) 



Y: m.*(A)(l-F(e;4j) = 1- E m,*iA)F{e\;^ 



AC[n] 



AC[n] 



^ 1 

n]\A> 



Proof. By ([2]), we have 



AQ[n\ ^ieA -I 

= ^M) Pr(ri >t\/ieA) 

= Y^M)Rief). 

AC[n] 

Similarly, using the dual Mobius form of 0^ (i.e., the fourth expression in Tabled]), we have 
Rs{t)=E[UX{t))] = 5] m,.(A)E[UX,(t) 

AC[n] '-ieA 

= Y (A) (l - Pr(Ti ^tyieA) 

AC[n] 

= 5] m,.(A)(l-F(e 

AQn] 



and for the last formula, we use the fact that J2Ac[n] ''ti„*(^) = 0t)*(l) = 1- 



□ 



It is noteworthy that Theorem [5] immediately provides concise expressions for the mean 
time-to-failure of the system, namely 



MTTF5 
MTTF 



Y m,{A) / R{e'/)dt 

AQn] ^0 

s = Y ^^'(A) / (l-F(ef; 



£,CXD 

n]\A 



dt. 



AC\n] 



Theorem [5] may suggest that the complete knowledge of the joint survival (or joint distri- 
bution) function is needed for the calculation of the system reliability function. Actually, as 
Theorem H] shows, all the needed information is encoded in the distribution of the indicator 
vector X(t). In turn, the distribution of X(t) can be easily expressed (see Dukhovny [2] and 
Dukhovny and Marichal ^) in terms of the joint probability generating function of X(t), 
which is defined by 



G{z,t) = E 



i=l 



{\zi\ ^l,iE [n]). 



As it is well known, the joint probability generating function has the advantage of being an 
expectation and yields not only the probabilities alone but also all kinds of moments via 
derivatives. 

By definition, we have 



G{z,t)= Y Pr(X(t)=x) n 4 

xe{0,l}" i&[n] 



Y Pr(X(t) = e'/) n 

ACfn.l ieA 



(14) 
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and hence G(z, t) is a multilinear polynomial in zi, Zn, which can be rewritten as 

G{z,t)= y: G{e'/,t)i[z, n (1-^^)- 

AC[n] i&A i£[n]\A 

Moreover, we can easily show [21 E] that G'(e^\f) = F{e^^). 
On the other hand, from f|T^ it follows that 

G{e'j^,t)= 5: Pr(X(t) = e^^) 

which shows that the set function A Pr(X(t) = e^^) is the Mobius transform of the set 
function A i— >• G(e^"^,t), that is, 

Pr(X(t) = e^^)= E(-l)'^'~"''G(e^\^)- 

B<ZA 

Combining this latter formula with ffTOj) enables us to express the system reliability function 
in terms of (^(z, t). 

3.2 Independent failures 

In the case when Ti,...,T„ are independent, which implies that the indicator variables 
. . . are independent for all t ^ 0, from (fT2|) we obtain the well-known formula 

Rs{t)= E v{A)Y[R,it) n il-R^it))■ (15) 

AC[n] iGA ie[n]\A 

Combining ([3]) and f|T5l) . we immediately retrieve the following classical formula (see for 
instance Rausand and H0yland [151 §4-5]) 

Rs{t)=^,{Ri{t),...,Rn{t)) 

and so both Rs{t) and MTTF5 can be expressed in different forms, according to the ex- 
pressions of (py chosen in Table [H For instance, using the primal Mobius form of (p^, we 
obtain 

AC[n] i&A 



MTTF 




3.3 Some examples 

Let us now examine some typical examples by considering both their structure functions and 
the corresponding l.p. functions: 

1. Series structure. If all the components are wired in series, we have 0^(x) = YliXi and 
P«)(t) = Aiti- In this case, w{A) = 6 if and only ii A = [n]. More generally, we can 
show [9l §5.4] that any l.p. function fulfilling the functional equation 

Pwih At[, . . . ,tn At'J = Pnj{ti, . . . ,tn) A Pw{t'i, ■■■,t'n) (17) 
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is of the form Pw(t) = AieB U for some subset -B C [n]. It then corresponds to a serially 
connected segment of components. 

The reliability of a series structure with n elements is given by Rs{t) = Pr(A"^^Tj > t). 
Using Theorems H] and \5[ we also have 

Rs{t) = Pr(X(t) = 1) = i?(ej;p = Pr(Ti > t, . . . > t). 

2. Parallel structure. If all the components are wired in parallel, we have 0t,(x) = Ui^i 
and Pwi't) = \Jiti. In this case, w{A) = 6 if and only if A ^ 0. Similarly to the series 
structures, we can show that any l.p. function fulfilling the functional equation 

V t;, . . . , t„ V 4) = . . . , t„) V pUt[, t'j (18) 

is of the form Pwit) = yieB^i some subset B C [n]. It then corresponds to a 
subsystem of parallel components. 

The reliability of a parallel structure with n elements is given by Rs{t) = Pr(V"^^Tj > 
t) . Using Theorems H] and [5l we also have 

Rs{t) = 1 - Pr(X(t) = 0) = 1 - F{e'^) = 1 - Pr(Ti ^ t, . . . ^ t). 

3. k-out-of-n structure, for some A; G [n]. By definition, a A;-out-of-?2 structure is charac- 
terized by the structure function 

1, if 'ZiXi ^ k, 
0, if J2iXi < k. 



It is then easy to show that 



u = n n^^ 

AC[n] i&A A<Z[n] ieA 



\A\=k |A|=n-fc+l 

which, in turn, entails 

Pw{t)= y /\ti = A y ti = /„-fe+i(t), 

AC[n] ieA A<Z[n] ieA 

\A\=k \A\=n-k+l 

where, for any A; G [n], : L" ^ L is the fcth order statistic function (see for instance 
Ovchinnikov [l2]). We recall [9], §5.5] that the n order statistic functions are exactly 
those l.p. functions that are symmetric in their variables. It follows immediately that 
a structure is of A;-out-of-n type for some k G [n] if and only if its system lifetime 
is a symmetric function (which is fn-k+i) of the component lifetimes. In this case, 
w{A) = 6 if and only if \A\ ^ k, which means that the system is functioning if at least 
k components are functioning. Clearly, the minimal representation ([S]) is such that 
Uw{A) = b if and only 1^41 = k. We also observe that 



m^. (A) 



0, otherwise, 

\_l)\A\-n+k-l^\A\^iy ii\A\^n~k + l, 

0, otherwise. 
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The reliability of a A;-out-of-n structure is given by Rs(t) = Pr(/„_fc_|_i(T) > t). Using 
Theorem m we also have 

n 

Rsit) = J2FTmt)\=j) = Pr(|X(t)| ^ k), (20) 

j=k 

and we can show [21 E] that Pr(|X(t)| = j) = [x^]G{xl,t) is the coefficient of in the 
nth degree polynomial G{xl,t). On the other hand, combining (fT3ll and (flQll gives 

j=k ^/ AC[n] 

\A\=j 



Example 6. When Ri{t) = e (z = l,...,n), it is convenient to calculate MTTF5 by 
using formula f|T6l) . Indeed, in that case, setting = Z^iga-^j) we simply obtain (see ^0\) 



MTTF5= E ^viA) / e-^^*dt= ^ 



^ ^ 21 



Assuming further that the structure is of A;-out-of-n type, by (HM we immediately obtain 



\A\-1\ 1 

AQn] \ ^ - 1 / 



MTTF5 = (-1)'^'"^ 

AC[n] 
\A\^k 



4 Systems with lower and upper bounds on lifetimes 

Consider now a more general system in which we allow upper and/or lower bounds on 
lifetimes of certain subsets of components. As shown by Dukhovny and Marichal [3], the 
structure of such a system can be modelled by means of a w.l.p. function, which is an l.p. 
function constructed from both variables and constants. 



4.1 System descriptions 

We first recall the definition of w.l.p. functions (see Goodstein |1] and Rudeanu [HI Chap- 
ter 3, §3]). 

Definition 7. The class of weighted lattice polynomial (w.l.p.) functions from L" to L is 
defined as follows: 

(i) For any k G [n] and any c & L, the projection t t^ and the constant function t 1— > c 
are w.l.p. functions from L" to L. 

(ii) If p and q are w.l.p. functions from L" to L, then p A q and p\/ q are w.l.p. functions 
from L" to L. 

(iii) Every w.l.p. function from L"^ to L is constructed by finitely many applications of the 
rules (i) and (ii). 
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It was proved [1] that any w.l.p. function p : L"' ^ L can be expressed in disjunctive and 
conjunctive normal forms, that is, there exist set functions w'^ : 2["'l — >■ L and w'^ : 2'"] L 
such that 

p{t)= V («;V)aA^O= a (w''{A)y\/tX 



AC[nl ^ igA ^ AC\n] ^ i&A 

.4 



Moreover, it can be shown |8] that, from among all the possible set functions w that 
disjunctively define p(t), only one is nondecreasing. Similarly, from among all the possible 
set functions w'^ that conjunctively define p(t), only one is nonincreasing. These special set 
functions are given by 

w\A)=p{ef) and w\A) = p{e;i^^^) . 

The w.l.p. function defined by a given nondecreasing set function w : 2'"] L will henceforth 
be denoted pw 

The following theorem, which generalizes Theorem [2] to w.l.p. functions, shows that the 
system is no longer characterized by a single structure function but by a one-parameter 
family of structure functions. 

Theorem 8. With any system fulfilling Ts = Pw{Ti, . . . ,T„) for some w.l.p. function p^ : 
^ L is associated a unique family of nondecreasing and nonconstant structure functions 
: t ^ 0}, where Vt{A) = lnd{w{A) > t) for all AC [n], such that 

Xs{t) = cP,,iMt)) (t^O). 
Proof. We follow the same reasoning as in the proof of Theorem [2l For any t ^ 0, we have 

Ind(p^(T) >t) = Ind( V (^(^) ^ f\T.i)>t 

AC[n] ieA 

= V (lnd(«;(A) >t)A/\ lnd{Ti > t) 

AC[n] ^ i&A 

= EvM)Ex,{t) 

AC[n] i£A 

= 0..(x(t)). 

Hence, we have Ts = p«,(T) if and only if Xs{t) = 0t,j(X(t)) for all t ^ 0, which completes 
the proof. □ 

Remark 9. According to Theorem [HI when modelling systems with collective bounds, it seems 
much more convenient to use w.l.p. functions rather than families of structure functions. 

The properties of the family of structure functions can be derived from the properties of 
the corresponding w.l.p. function. Let us examine some of them: 

1. Boundary conditions. We have 0„t(O) = Vt(0) = lnd{w{0) > t) and <pvt{'^) = 
Vt{[n]) = lnd{uj{[n]) > t). 

2. Pivotal decomposition. The median-based decomposition formula ([7]), which also 
holds for any w.l.p. function, leads again to the pivotal decomposition of each structure 
function 6,,/. 



X) = Xi(j)y^{li,X) + (1 - Xi)(Py,{Oi,X) (t ^ 0). 
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3. Minimal representations. From a nondecreasing set function w : 2["1 L, define 
the set functions ut : 2'"] — > L and -uf,, : 2["'l — L as 



w{A), if w{B) < w{A) for all B ^ A, 

a, otherwise, 

'w{A), if w{A) < w{B) for all B ^ A, 

b, otherwise. 



The disjunctive and conjunctive representations of the w.l.p. function having a 
minimal number of terms write (see Marichal P, Proposition 8]) 

pM= V (ui{A)Af\u)= A (<{[n]\A)v\/t, 

AC[n] ^ ieA ^ AC[n] ^ iGA 

The corresponding form of the family of structure functions follows. 

4. Dual structure. The dual of a family \^(^vt '■ t ^ 0} of structure functions is the 
family {0^,* : t ^ 0} of structure functions such that 

vl{A) = 1 - lnd{w{[n] \A)>t). 

It follows that there is no set function w* : 21"] L such that v^{A) = Ind(«;*(A) > t). 

Example 10. Consider the bridge structure as indicated in Figure [1] and assume that the 
time-to- failure of the central component must lie in the interval [I, u] for some ^ I ^ u ^ oo. 
We immediately see that the minimal path sets are Pi = {1,4}, P2 = {2,5}, P3 = {1,3,5}, 
and P4 = {2,3,4}. Hence, the w.l.p. function associated with this structure is given by 

Pu)it) = {ti A ^4) V (^2 A ts) V (ti A median(/, t^, u) A t^) V (t2 A median(/, ^3, u) A ^4) 

= {ti A U) V (ta A tg) V (/ A ti A U) V (m A ti A ta A U) V (/ A A U) V (w A A A ^4) 

and the corresponding family of structure functions is 

0^((x) = (X1X4) n (X2X5) n (Ind(/ > t)a:iX5) n (Ind(M > i(:)xiX3X5) 
n (Ind(/ > t)2;2a;4) 11 (Ind(u > t)x2X^X4). 




Figure 1: Bridge structure 
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4.2 Exact reliability formulas 

Regarding the reliability calculation, Dukhovny and Marichal [3] established the following 
result, which is a direct generalization of Theorem |H 

Theorem 11. We have 

Rs{t) = E vM)'Pr{X{t) = ef), (22) 

AC[n] 

Mt) = 1- E <(A)Pr(X(t) = e|'4j. 

AC[n] 

Similarly, a direct generalization of Theorem [5] is stated in the following theorem: 
Theorem 12. We have 

Mt) = E rn,M)RieA), (23) 

AC[n] 

Rsit) = E (l - n%]\A)) = 1 - E (24) 

AC[n] AC[n] 

As far as the mean time-to-failure of the system is concerned, from fl23l) and fl24l) we 
immediately obtain 



rw(B) 

MTTF5 = E E(-l)'^'""" / Rie'/)dt, (25) 
MTTF, = E E(-l)'^''"'' rn«.(^-^(%M))^^- (26) 



AC[n] -BCA 



When the variables Ti, . . . ,r„ are independent, from fl22p we immediately retrieve the 
formula (see Marichal [11]): 

Rsit)= E v,iA)l[R,{t) n (l-^.(^))- 

AC[n] ieA ie[n]\A 

Considering the family {(p^^ : t ^ 0}, where is the multilinear extension of (p^^, we then 
observe that 

Rs{t)=^jR^{t),...,RrXt)) (27) 

and (f)^^ can be chosen from among the forms given in Table [H where each v should be 
replaced with Vt- Also, from fl25p and (1261) we immediately derive 

rw(B) 



MTTF^ = E E(-l)'^'""" / \{R:{t)dt, (28) 

(•OO 

MTTYs = E E(-l)'^'""'' / n^.(t)rft. 

AC[n] BCA ^«)([n]\S) ieA 

Let us now examine some examples by considering special w.l.p. functions. They gener- 
alize the classical examples considered in §3.31 (series, parallel, and fc-out-of-n structures). 
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1. Weighted minimum. A weighted minimum function is a w.l.p. function '■ L"^ ^ L 
whose underlying set function t/; : 2["1 — L fulfills 

w{[n] \ (A U B)) = w{[n\ \ A) A w{[n] \ B). 

Such a function fulfills equation ([T7|) and is of the form (see P, §5.2]) 



■n 

Pu,{t) = f\[w{{i])yu). 



i=l 



It then corresponds to a series structure with a lower bound on the lifetime of each 
component. By using (12^ . we can easily show that 

Rs{t) = i- E (-i)''"^'^(ef„r\A)n(i-^*(w)) 

AC[n] i€A 

and, in case of independence (see f l27|l ). 

n 

Rsit) = l[{vt{{t})UR,it) 

i=l 



2. Weighted maximum. A weighted maximum function is a w.l.p. function p^; : L"^ ^ L 
whose underlying set function w : 2^'^'^ ^ L fulfills 

w{AuB) = w{A) yw{B). 

Such a function fulfills equation f[T5]) and is of the form (see P, §5.2]) 



7t 

pM = V («^(W) Ati). 



i=l 



It then corresponds to a parallel structure with an upper bound on the lifetime of each 
component. By using (123!) . it is also straightforward to show that 

AC[n] i&A 
A^0 

and, in case of independence (see f l27|) ). 

n 

Rs{t) = l[^ti{^})R^it). 

i=l 



3. Symmetric w.l.p. function. We can generalize the k-out-oi-n type structures simply 
by considering symmetric w.l.p. functions p^ '■ L"^ ^ L. The underlying set functions 
are cardinality based, i.e., such that w{A) = w{B) whenever \A\ = \B\. If we define 
the function w : {0, 1, . . . , n} — L by w{A) = w{\A\), we can easily show [3] that any 
symmetric w.l.p. function can always be put in the form 

n 

Pw{t) = V (*(^) A /n-fe+l(t)), 
k=0 
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where fn-k+i is the order statistic function defining the fc-out-of-n structure (see §3.31) . 
Moreover, we can show that 

Rs{t) = Pr(/„_,(,)+i(T) >t)= Pr(|X(t)| ^ k{t)), 

where k{t) := min{A;,n + 1 : w{k) > t}, which generahzes f l2UI) . For more details, see 
Dukhovny and Marichal [3]. 

Example 13. Let us calculate MTTF5 when Ri{t) = e~^»* {i = l,...,n). Using ([28]) and 
setting Xa = J^ieA^iy we obtain (see [TT] ) 



AQ[n] BCA 

E E(-i)'^'""'' r- 

AC[n] BCA 

Ayt0 



rw(B) 

MTTFs = E E(-l)'^''""/ ^'^""'dt 

I _ g-AAw(-B) 



When Pu, is an l.p. function (that is, w{A) G {0, 00} and w{0) = 0), we retrieve formula 
(EH). Indeed, 

MTTF5= E l^ E (-1)1-^1-1^1= E 

^7^0 w{B)=oo A^0 



5 Distribution functions of w.l.p. functions 

The articles [2|, [TOl [11] on which this paper is partly based were motivated by the exact com- 
putation of the distribution functions and the moments of l.p. functions and w.l.p. functions 
of random variables. 

In this final section we point out the fact that calculating the distribution function of 
an arbitrary w.l.p. function amounts to calculating the reliability function of a semicoherent 
system with possible lower and upper bounds on component lifetimes. 

Let L C M be a totally ordered lattice, let : L" — > L be a w.l.p., and let Ti, . . . , T„ be 
L-valued random variables. 

The distribution function of the random variable p«,(T) is defined as 

Fpjt) = Pr(p^(T) ^ t) (teL). 

Clearly, this function fulfills the identity Fp^{t) = 1 — Rs{t), where Rs{t) is the reliability 
function of the coherent system described by the w.l.p. function p^,. 

Using formulas (l22|) - f[M[) . we then obtain immediately the following formulas for Fp^[t): 





it) 


= 1 - E PKX(t) = e^^) 

AC[n] 


Fv 


it) 


= E <(^)Pr(X(t)=ej'4j 

AC[n] 


Fpw 


it) 


= 1- E m^AA)R{e/) 

AC[n] 


Fpyj 


it) 





AC[n] 
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where Vt{A) = lnd{w{A) > t) and v*{A) = 1 - lnd{w{[n] \A) > t). 

When the arguments Ti,...,r„ are independent, each Tj having distribution function 
Fi{t), we obtain (see Marichal 



1- E MA) m^-m) n m) 



AC[n] i&A ie[n\\A 



AC[n] ieA ie[n]\A 



FpSt) 



1- ^^M) m^-F-m 



FpSt) 



AC[n] iGA 



6 Conclusion 

We have discussed the formal parallehsm between two representations of semicoherent sys- 
tems: structure functions and Lp. functions. Their languages are shown to be equivalent in 
many ways. The l.p. language is demonstrated to have significant advantages. One is the 
natural generalization to w.l.p. functions and corresponding systems with bounded subsys- 
tem lifetimes. The other is the fact that, due to the distributive property of the indicator 
function Ind(-) with respect to lattice operations (see proofs of Theorems [2] and [8]) , the l.p. 
description is a very natural tool to connect the system's structure to the lattice of typical 
reliability events of the kind T ^ t, to connect the system's purpose, as encoded in the l.p. 
function, to the system's equipment, as expressed in the joint distribution of units' lifetimes. 
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